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PREFACE 

Work  on  the  development  of  a  long  base-line  interferometer  for 

wavelengths  in  the  neighborhood  of  1  cm.  started  in  the  Space  Sciences 

* 

and  Radioastronomy  Laboratories  in  1970.  In  our  annual  report  for 
1970-1971  we  described  the  interferometer  which  we  had  built  at  the 
Hat  Creek  Radioastronoray  Observatory  and  noted  that  radio  sources  had 
been  located  which  could  serve  for  highly  accurate  work  on  geodesy 
and  navigation.  Prof.  Welch  (co-principal-inves tigator)  reported  the 
work,  early  in  1972,  to  the  annual  TAC  meeting  for  the  review  of  the 
electronics  program  held  by  the  Electronics  Research  Laboratory  of 
the  College  of  Engineering  and  at  that  time  put  forth  our  ideas  about 
the  use  of  an  interferometer  on  board  a  ship  as  an  instrument  for 
navigation. 

We  then  proceeded  to  study  this  matter  seriously.  The  following 
report  covers  a  rather  comprehensive  study  carried  on  by  Mr.  Michael 
Morgan,  a  graduate  student  in  the  Dept,  of  Electrical  Engineering  and 
Computer  Science,  under  the  direction  of  Prof.  Silver.  The  report  is 
not  only  a  feasibility-study  but  also  delineates  the  techniques  to 
be  employed,  both  analytically  and  operationally,  for  locating  the 
position  of  a  ship  on  the  sea  from  a  knowledge  of  the  position  of  a  radio- 
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source  in  the  sky. 

For  the  wavelength  region  of  1  cm.  a  long  base-line  interferometer 
can  be  established  on  board  a  ship  by  placing  the  two  elements  —  an¬ 
tennas  —  at  the  fore  and  aft  positions  of  the  deck;  all  of  the  lines 
and  complex  electronics  can  be  located  below  deck;  it  is,  in  fact, 
desirable  to  do  that.  It  is  possible  to  go  a  step  further,  for,  un¬ 
less  it  is  necessary  to  chart  position  continuously,  the  antennas  can 
be  used  for  other  operations.  One  needs  only  a  switch  arrangement 
whereby  the  antennas  can  be  connected  to  operate  in  the  interferometer¬ 
mode  or  be  disengaged  from  the  electronics  of  the  interferometer  and 
tied  in  to  other  systems.  The  point  of  mentioning  this  is  to  allay 
arguments  about  costs  and  that  the  interferometer  will  rob  the  9hip 
of  other  capabilities. 

Whether  or  not  the  techniques  developed  in  this  report  find  their 
way  into  naval  operations  depends  on  the  evaluation  made  of  them  by  the 
Department  of  the  Navy.  There  are  undoubtedly  many  operational  consider¬ 
ations  to  be  taken  into  account  but  such  matters  lie  outside  our  area 
of  competence  and  responsibility.  The  important  point  is  that  this  work 
stands  as  an  example,  of  considerable  merit  in  this  writer's  opinion, 
of  how  the  pursuit  of  basic  research,  in  this  case  in  radioastronomy, 
can  and  does  lead  to  applications  and  operational  systems.  We  trust 
that  what  we  are  presenting  in  this  work  does  have  practical  utility. 

'  Samuel  Silver 

'  Principal  Investigator 
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/  I.  INTRODUCTION 

/ 

The  subject  of  this  report  is  an  examination  of  the  feasibility 
of  using  a  long-baseline  interferometer  aboard  a  ship  as  a  navigational 
device.  The  ideal  theoretical  resolution  of  such-  an  interferometer  with 
a  correlation-receiver  has  the  potential  of  determining  the  latitude 
and  longitude  of  a  ship  at  sea  with  great  accuracy.^  For  example,  an 
interferometer  whose  base-line  is  200  meters,  operating  at  a  wavelength 

of  1  cm.,  has  a  central  fringe  of  about  10  seconds  of  arc.  By  observing 

i 

discrete  (i.e.,  point-like)  radio  stars,  whose  celestial  coordinates 
have  been  determined  by  land-based  observatories,  in  the  near  zenith 
position  relative  to  the  ship  the  position  of  the  ship  on  the  earth  can 
be  established  with  an  angular  accuracy  of  about  the  resolving  power  of 
the  interferometer.  For  the  example  quoted, a  conservative  estimate  of 
this  accuracy  is  10  sec.  of  arc  for  the  angular  coordinates  which  re¬ 
presents  location  of  position  to  an  accuracy  better  than  1000  feet. 

In  addition  to  the  advantage  of  resolving  power  there  is  the  advan¬ 
tage  that  the  instrument  is  not  limited  to  nighttime  use.  As  in  the 
optical  case,  however,  the  use  of  the  system  Is  limited  by  meteorological 
conditions.  The  limitations  imposed  by  meteorological  conditions  are 
strongly  dependent  on  the  frequency  (or  wavelength)  at  which  the  system 
operates . 

The  report  that  follows  has  two  major  subdivisions.  The  one 
deals  with  topics  —  Preliminary  Topics  —  that  provide^  the  framework 

for  computing  position  from  certain  types  of  data;  the  >other  deals  with 

A  ' 

the  special  considerations  of  interferometry  pertinent  to  a  navigational 
system  and  the  navigational  technique  proper. 

•'/X 


II.  PRELIMINARY  TOPICS 


2.1  Theory  and  Measurement  of  Ship  Motions 

The  navigational  system  envisaged  in  this  report  consists  of  a 
microwave  interferometer  mounted  on  board  a  ship.  On  the  ocean,  the 
ship  is  a  moving  platform  whose  motions  are  made  up  of  translation 
and  rotations.  The  effect  of  such  a  moving  platform  on  the  operation 
of  the  interferometer  is  one  of  the  major  considerations  in  the  devel¬ 
opment  of  the  navigational  system.  For  this  reason  we  present,  in 
this  section,  some  elementary  information  from  the  theory  of  ship 
motions.  A  much  more  detailed  development  can  be  found  in  such  ref¬ 
erences  as  those  by  Blagoreshchenskii,  [1],  or  Beurs,  [2]. 

A  ship  on  the  ocean  is  sketched  in  Figure  1.  The  cartesian  axes 
x  ,  x  ,  x  are  the  body-fixed  principal  axes  of  the  ship.  The  general 

123  i. 

motion  of  the  vessel  can  be  resolved  into  translation  of  the  origin 
which  is  assumed  to  be  at  the  center  of  gravity  of  the  vessel  and  ro¬ 
tations  about  the  principal  axes.  The  angular  displacements  are  shown 

in  the  figure  as  £  (t),  f;  (t) ,  5  (t);  in  nautical  circles  they  are 

1^3 

known  as  the  roll,  pitch,  and  yaw  angles,  respectively.  It  is  known 
from  experiment,  [1;  pp.  338-351],  that  wave  height  on  the  ocean  and, 
consequently,  the  roll,  pitch  and  yaw  angles  have  amplitudes  that  are 
random  variables  having  Gaussian  distributions.  'Fhe  typical  sample 
function  for  any  of  the  three  random  variables  £^(t),  i  =  1,2,3,  is 
sketched  in  Figure  2.  The  time  dependence  of  these  rotational  dis¬ 
placements  in  mild  and  moderate  sea-states  is  quasi-sinusoidal  with 
an  almost  constant  period  of  motion  and  a  Gaussian  distribution  of 


amplitudes.  Typical  standard  deviations  of  roll,  pitch,  and  yaw 
angles  as  functions  of  seastate  for  an  800  ft.  vessel  are  given  in 
Table  1.  Some  typical  periods  of  rolling  in  a  mild  seastate  for 
various  size  vessels  are  given  in  Table  2. 

The  geometrical  factors  relating  the  interferometer  and  a  radio 
source  to  the  principal  axes  are  shown  in  Figure  3.  The  base-line 
of  the  interferometer  is  along  the  x  -axis  and  its  direction  is  speci¬ 
fied  by  the  unit  vector  3(t)  which  in  a  space-fixed  reference  system 
is  a  function  of  time.  The  radio-star-source  is  located  by  the  unit 
vector  t(t)  which  is  a  function  of  time  in  both  the  ship-fixed  axes 
and  the  space-fixed  axes  by  virtue  of  the  motion  of  the  ship  and  the 
rotation  of  the  earth.  The  angle  between  s(t)  and  r(t)  is  denoted  by 
0(t),  The  length  of  the  interferometer  baseline  is  given  by  S(t). 

There  will  be  time-dependent  perturbations  in  this  baseline  length  due 
to  both  the  effects  of  mechanical  stress  and  thermal  expansion.  To 
obtain  the  maximum  achievable  resolution  of  the  interferometer,  of  the 
order  of  one  "fringe  width"  it  is  necessary  to  know  the  length  of  the 
physical  baseline,  at  all  times  during  the  operation  of  the  system, 
to  within  a  fraction  of  a  wavelength.  The  variation  in  length  due  to 
mechanical  stress  will  depend  upon  the  construction  of  the  vessel.  The 
order  of  magnitude  of  thermal  effects  can  be  obtained  by  assuming  that 
the  ship  is  constructed  of  steel  with  a  mean  linear  coefficient  of  ex¬ 
pansion  of  1.2  x  10-5  (C°)-1,  A  temperature  variation  of  25  °C  =  45  °F 
on  a  vessel  made  of  such  material  generates  a  change  in  length  in  the 
amount  of  6  cm.  in  a  baseline  whose  reference-length  is  200  meters.  If 
the  wavelength  of  operation  is  1  cm. ,  the  apparent  location  of  sources 


Figure 


Figure  2 


Ship  Principal  Axes  and  Ro 
Angles  v 


Typical  Sample  Function 


A 


A 


A 


Displacement 


,*2 

rc2(t) 


A--  -  df.  standard 
deviation 


7 


Table  1  Typical  Standard  Deviations  (800  ft. 


Seastate 


0 

2 

3 

6 

Roll  a? 

2 

0° 

1-2° 

5-6° 

>10 

Pitch  o£ 

2 

0° 

2-3° 

3-4° 

4-6° 

Yaw  aE, 

0° 

0-1° 

3-4° 

5-8° 

3 


vessel) 


Table  2 _ Typical  Periods  of  Rolling  (mild  seastate) 


Vessel  Type 

Displacement  (Tons) 

Period 

Ocean  Liners 

30,000-50,000 

20-28 

Battleships 

<  10,000 

14-18 

Heavy  Cruisers 

<  10,000 

13-15 

Destroyers 

<  10,000 

7-9 

Patrol  Boats 

<  10,000 

4-6 
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can  be  changed  by  as  much  as  six  fringe  widths.  It  is,  therefore, 
necessary  to  have  a  subsystem  for  monitoring  accurately  the  length  of 
the  baseline.  Several  techniques  are  available  for  this  purpose;  one, 
for  example,  is  to  measure  pulse  return  time  on  a  taut,  but  stretchable, 
shorted  transmission  line.  We  will  assume  in  the  subsequent  discussion 
that  S(t)  is  known  to  within  a  fraction  of  a  wavelength  at  all  times. 

The  technique  for  determining  the  roll,  pitch,  and  yaw  angles  of 
the  interferometer  baseline  utilizes  highly  accurate  accelerometers, 
placed  judiciously  on  the  vessel,  as  shown  in  Figure  4.  The  sensitivity- 

l 

directions  of  these  linear  accelerometers  are  shown  by  the  attached 
vectors.  The  magnitude  of  linear  acceleration  is  given  by  a^(t) 
or  b^(t),  i  =  1,2,3.  The  subscript  (i)  is  related  to  the  angular  dis¬ 
placement  derived  from  the  particular  accelerometer’s  data.  The 

first  set  of  accelerometers  is  located  along  the  x  axis,  one  on  the 

2 

starboard  and  one  on  the  port  side  of  the  ship.  The  second  and  third 

sets  of  accelerometers,  (a  ,b  )  and  (a  ,b  ),  are  located  at  the  antenna 

2  2  3  3 

positions.  The  roll,  pitch  and  yaw  angles  are  obtained  by  "integrating" 
the  differential  equations 


S(t) 


2 


ds 

dt 


(1) 


for  t  =  2,3  with  appropriate  initial  conditions  on  (t).  The  equation 


for  £  (t)  contains  the  distance  W(t)  between  the  accelerometers  a  and 
1  1 

dW  dgi 


b  and  there  is  no  cross-product  term  _ 

1  at 


dt 


The  extent  to  which  the  time  dependence  of  the  length  of  the  base¬ 
line,  S(t),  and  the  time  dependence  of  W  is  significant  depends  on  thermal 
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and  mechanical  stresses.  Data  on  these  factors  will  have  to  be  ob¬ 
tained  before  we  pass  judgment  on  the  question. 

To  a  first  order  approximation  in  which  we  neglect  the  cross¬ 
term  ^  -^T  the  angles  £^(t)  are  given  by 


n> .  f y  » 

1  Jo  Jo  H  <") 


(n)  3 


dq  dx 


(2) 


e  (t) 


a 


(a2(n)  -  b2  (n) ] 

s  (n) 


dn  dx 


(3) 


e  (t)  = 


■u 


[a^Cn)  -  b3(n) ] 
S  (n) 


dn  dx 


(4) 


on  the  assumption  that 


q<o) 


(0)  =  0 


i  -  1,2,3 


(5) 


The  desired  accuracy  of  the  measured  roll  angle,  £^(t),  is  not 
extremely  critical  because  rotation  about  the  interferometer's  base¬ 
line  does  not  affect  the  fringe  pattern.  Rolling  motion  will,  however, 
shift  the  position  of  the  source  in  the  main  beams  of  the  element- 
antennas  of  the  interferometer.  The  angles  and  are  needed 
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to  slowly  track  the  two  antennas  to  compensate  for  the  ship's  motion 
and  keep  the  position  of  the  source  fairly  constant  in  their  main  beams. 

For  this  reason  the  error  in  the  determination  of  £  j(t)  should  be  much 
smaller  than  the  main  beam-width.  Typical  beam-widths  at  centimeter 
wavelengths  for  moderate  size  parabolic  reflectors  are  usually  of  the 
order  of  a  few  minutes  of  arc.  For  a  5  meter  reflector  and  X  *  1  cm 
we  have  0^  7  min.  The  required  accuracy  of  the  pitch  and  yaw 

angles  is  much  more  stringent.  These  angles  are  used  not  only  in  the 
antenna  tracking  scheme  but  also  in  the  data-processing  scheme  for  deter¬ 
mining  tbe  ship's  position  from  the  measurements  of  the  phase-correlation 
of  the  waves  from  the  source  detected  by  the  two  antennas  of  the  inter¬ 
ferometer.  The  error-tolerance  of  the  pitch  and  yaw  angles  must  be  less 
than  a  fringe  width,  which  at  centimeter  wavelengths  is  of  the  order  of 
magnitude  of  several  seconds  of  arc.  The  use  of  highly  accurate  accelero¬ 
meters  located  at  the  antennas  should  give  £2(t)  C3(t)  with  very 

small  error.  Note  that  in  this  procedure  that  we  have  just  sketched 
even  the  effects  of  bending  of  the  baseline  due  to  lack  of  perfect  rigidity 
of  the  vessel  are  included  in  the  measurements. 

2.2  Coordinate  Systems  and  Transformations 

We  need  to  use  three  fundamental  coordinate  systems  in  developing 
the  theory  of  navigation  by  means  of  the  interferometer.  These  systems 
which  we  name  (a)  the  local  horizon  system,  (b)  the  latitude-longitude 
system  and  (c)  the  celestial  system  are  shown  in  Figures  5,  6,  and  7, 
respectively.  In  the  local  system,  Figure  5,  the  x' ,  y',  and  z' 

cartesian  axes  are  oriented  in  the  due  West,  due  South,  and  zenith  direc- 
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tions ,  respectively.  The  orientation  angles  of  the  unit  baseline 
vector,  s’ (t) ,  are  given  by  the  azimuth  and'  elevation  angles,  A(t) 
and  E(t),  respectively,  as  shown  in  Figure  5.  These  angles  should 

'  j  a 

be  monitored  at  all  times  and  are  given  by 

A(t)  =  5,  (t)  +  it  +  p  (t) 

0  O 

E(t)  .  =  (t) 

where 

(t)  =  pitch  angle 

S3  (t)  =  yaw  angle 

pQ  (t)  =  ship  bearing  (East  of  North) 

The  "apparent"  source  position  unit  vector,  ?a'  (t) ,  , points  in  the 
direction  that  the  source  appears  to  be  located  with  respect  to  the 
ship.  This  apparent  source  position  may  differ  from  the  actual  source 
position,  represented  by  the  unit  vector  £■'(!),  because  of  atmospheric 
effects.  These  effects  will  be  considered  in  more  detail  in  the  fol¬ 
lowing  section  where  an  approximate  linear  transformation  between 
t'(t)  and  f,a'(t)  will  be  derived  by  means  of  ray-theory  onv  the  basis 
that  the  effect  is  ascribable  to  refraction  by  the  atmosphere.  Note 
that  the  exact  cartesian  components  of  ?'(t)  or  £a' (t)  are  not  known 
a  priori  in  the  local  horizon  system. 


(6> 

(7) 
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In  the  latitude-longitude  coordinate, system.  Figure  6,  the  xy 

plane  is  the  equatorial  plane  of  the  earth  while  the  yz  plane  is  the 

l 

Greenwich  meridian  plane  and  the  z  axis  is  the  polar  axis.  The 
orientation  angles,  5  and  A,  of  the  unit  vector,  z',  in  the  direction 
of  the  zenith  relative  to  the  ship,  are  the  latitude  and  longitude 
angles,  respectively,  of  the  position  of  the  ship  on  the  earth.  It 
is  these  angles  which  the  interferometer  navigational  system  will 
supply  to  a  very  high  degree  of  accuracy.  The  unit  vector,  s^(t), 
is  parallel  to  the  x’y'  plane  and  corresponds  to  the  rest  position 
that  the  unit  baseline  vector,  S'(t)',  would  have  on  "a  perfectly  calm 
sea.  The  azimuth  and  elevation  angles  of  this  unit  baseline  "rest 
position"  vector  are  given  by 


- 

Vc) 

TT  +  PQ(t) 

(8) 

Eo(t) 

o  ; 

(9) 

The  cartesian  components 

of  §'(t)  and  §’ (t)  will  be  known  at  all 
o 

times 

during  the 

operation 

of 

the  interferometer-navigational-system  and  are 

given  by 

Sx(t) 

a 

sin  A(t)  cos  E(t) 

(10) 

s;(t) 

- 

cos  A(t)  cos  E(t) 

(ID 

sz(t> 

= 

sin  E(t) 

(12) 
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and 


s'  (t)  =  sinA(t) 

xo  ° 

Sy  (t)  =  COS  A^(t) 

s'  (t)  =  0 

2o 


(13) 

(14) 

(15) 


The  equivalent  latitude  and  longitude  angles  of  the  actual  source 
position  are  given  by  and  Agl  respectively.  These  angles  can  be 
found  at  any  time  using  the  tabulated  declination  and  right  ascension 
of  the  particular  source  being  considered.  The  cartesian  components 
of  the  unit  position  vector  of  the  actual  source,  f(t),  can  be  obtained 
using 


r  (t)  =  sin  A  (t)  cos  & 

x  s  s 


(16) 


r  (t)  =  cos  A  (t)  cos  6 

y  s  s 


(17) 


r  (t)  =  sin  6 

z  s 


(18) 


Note  that  the  exact  cartesian  components  of  the  unit  baseline  vectors 
3(t)  and  3Q(t)  are  not  known  a  priori.  It  will  be  shown  that  simulta¬ 
neous  knowledge  of  the  cartesian  components  of  any  vector  in  both  the 
local  horizon  system  and  the  latitude-longitude  system  is  equivalent 
to  knowledge  of  the  position  of  the  ship  on  the  earth. 


The  celestial  coordinate  system  is  shown  in  Figure  7.  The 
x"-axis  points  in  the  direction  of  the  vernal  equinox  which  is  also 
referred  to  as  the  "first  point  of  Aries."  (Years  ago  this  axis  was 
oriented  towards  the  constellation  Aries  while  at  present,  due  to 
precession  of  the  earth's  polar  axis,  the  x"  axis  points  towards  the 
constellation  Pisces.)  The  x"y"  plane  is  the  equatorial  plane  and  the, 
z"  axis  is  colinear  with  the  z  axis  of  the  latitude-longitude  system. 
The  position  angles  of  the  actual  source  location,  a  and  5  ,  are  the 
right  ascension  and  declination  angles,  respectively.  The  angles  will 
be  known  for  any  particular  source  using  tabulated  data  obtainable 
from  ground-based  measurements  of  positions  and  appropriate  correction 
procedures  for  the  given  epoch.  A  detailed  description  of  the  various 
standard  coordinate  systems  and  the  procedure  for  making  the  necessary 
corrections  in  the  position  of  a  source  for  a  given  epoch  using  the 
ephemeris  are  given  in  many  texts  on  optical  or  radio  astronomy  such  as 
that  by  Kraus,  [3,  chapter  2], 

The  declination  angle,  6  ,  and  the  equivalent  latitude  angle,  5  , 

s  s 

are, of  course, the  same.  The  relationship  between  the  right  ascension 

angle,  a  (t),  and  the  equivalent  longitude  angle,  X  (t) ,  is  shown  in 
s  s 

Figure  8.  The  angle,  y(t),  measured  in  radians,  is  given  by 

U(t)  -  Y2  (S-T‘  +  E’T‘)  (19) 

where 


S.T.  =  sidereal  time 


E.T . 


ephemeris  time 
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S.T.  can  be  obtained  directly  from  "The  American  Ephemeris  and  Nautical 
Almanac"  or  similar  ephemeris.  E.T.  can  be  obtained  from  Greenwich  mean 
time,  G.M.T.,  using 


E.T.  =  G.M.T.  +  AT  (20) 

where  the  correction  term,  AT,  is  given  for  the  particular  epoch  in 

any  standard  ephemeris.  The  relationship  between  a  and  X  (t)  is, 

s  s 

as  shown  in  Figure  8, 

a  +  X  (t)  =  y(t).  (21) 

s  s 

Later  in  this  report  we  will  use  explicitly  the  linear  transfor¬ 
mation  between  the  cartesian  components  of  a  vector  in  the  latitude- 
longitude  system  and  those  of  the  same  vector  in  the  local  horizon 
coordinate  system.  This  transformation  can  be  represented  in  matrix 
from  by 

U '  =  A  •  U  (22) 

where  U'  and  U  are  column  vectors  whose  elements  are  the  components  of  the 
given  vector  in  the  respective  coordinate  systems.  The  elements  of  the 
transformation  matrix  can  be  obtained  easily  by  considering  the  total 
transformation  to  be  composed  of  two  successive  unitary  rotation  trans¬ 
formations.  Using  the  1-2-3  representation,  x  ■«-»■  1 ,  y  ■«-*■  2 ,  z  3, 
the  transformation  in  (23)  can  be  written  as 


gure  7  Celestial  Coordinates 


Actual 

N.Pole  ^Source 

n  / 


gure  8  The  Equatorial  Plane 
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3 

u'  =  L  A  U  ,  i  =  1,  2,  3  (23) 

j  =  l  3  3 


The  elements  of  the  transformation  matrix.  A,  are  given  directly  in 
terms  of  the  latitude-  and  longitude-coordinates  of  the  ship,  6 ( t ) 
and  A(t),  respectively. 


A(t)  = 


cos  A(t) 


-sin  A(t) 


sin  A(t)  sin  <5(t)  cos  X  ( t )  sin  <5(t)  -cos  6(t) 


Lsin  A(t)  cos  6 ( t )  cos  X(t)  cos  6 ( t )  sin  6  (t)J 


(24) 


As  was  previously  mentioned,  the  simultaneous  knowledge  of  the  car¬ 
tesian  components  of  any  vector  in  both  the  horizon  system  and  lati¬ 
tude-longitude  system  is  equivalent  to  a  knowledge  of  the  ship's  lo¬ 
cation  on  the  earth.  We  will  now  obtain  equations  of  evolution  of 
the  terms  in  this  transformation  matrix.  These  equations  will  be 
needed  at  a  later  time.  We  note  that  the  latitude  and  longitude  angles 
may  be  expressed  in  the  form 


5(t)  =  6 (tQ)  +  A5  (t) 


(25) 


A (t)  =  A ( t  Q )  +  AA(t) 


(26) 


The  changes  of  these  angles  may  be  obtained  by  considering  the  motion 
of  the  vessel  on  the  planar  projected  latitude-longitude  curvilinear 
coordinates  shown  in  Figure  9.  The  equations  of  evolution  of  these 
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position  coordinates  of , the  ship  are  given  by 


tf 


<5(t)  =  5  (tQ)  +  -  I  v(t)  cos  pD  (t)  dr 


(27) 


X(t)  = 


Ht0) 


-il 


v(t)  sec  6(t)  sin  p0(x)  dx 


(28)- 


where  pQ(t)  =  measured  bearing  of  ship  (E  of  N) ,  (radians) 


v(t)  =  measured  speed  of  the  ship,  (m/sec) 


radius  of  the  earth,  (meters) 


We  now  substitute  (27)  and  (28)  into  the  matrix  elements  of  (24) 
and  use  the  angle  sum  formulas:  Defining 

t 


f.tt)  £  cos  I 


/ 


v(t)  sec  6(t)  sin  pQ(x)  dT 


(29) 


f2(t)  =  sin 


in  (»  / 


v(t)  sec  6 ( t )  sin  pQ(x)  dx 


(30) 


and 


g^t)  =  cos 


f 

il 


v(x)  cos  p0(x)  dx 


(31) 
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^  f 

iR  J 


g2(t)  =  sin  —  J  v( t )  cos  p0(i)  dr  j 


(32) 


and  using  the  abbreviations 


a.  .  (t  )  =  a  ,  . 

iJ  0  1J0 


i,  j  =  1,  2,  3 


(33) 


we  obtain  the  equations  of  evolution 


31  1  (t)  =  fiail  +  f2a!2 


(34) 


ai2(t)  fiai2o  “  f2ailo 


(35) 


a21(t)  =  f igia210  ‘  flg2a31Q  +  f2gia220  “  f282a320 


(36) 


l22(t)  flSla22o  f 1 ®2a 3 2Q  ~  ^2^1 a2 1 0  +  f2S2a31Q 


(37) 


a  (t:)  gla23  82a33 
3  o  o 


(38) 


a31(t).  =  flgia310+fl82a21n+f2gia32n+f2g2a220 


(39) 


a32(t)  =  flgia32  +  flg2a220  "  f2gia3l0  '  ^2*21 

o  o  o  o 


(40) 


8la33  +  g2a23 

o  o 


a33(t) 


(41) 
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2.3  The  Refraction  Matrix  Approximation 

As  was  mentioned  in  the  previous  section  we  will1  now  derive  an 
approximate  linear  matrix  transformation  between  the  unit  vector,  , 
r'  (t),  in  the  direction  of  the  actual  source-position  and  the  unit 
vector,  r  ' (t)  ,  in  the  direction  of  the  apparent  source-position, 
using  ray  theory.  The  ray  path  in  the  earth's  atmosphere  of  an  inci¬ 
dent  monochromatic  plane  wave  of  frequency  id  is  shown  in  Figure  10. 

The  curved  ray  path  taken  by  the  incident  field  will  be  determined 
by  Fermat's  principle  [4,  p.  356]. 

6  J  n(r,  id)  ds  =  0  (42) 

s 

where 

n(?,  oj)  =  index  of  refraction  as  a  function  of  position 
and  frequency. 

If  n(r,  m)  is  a  known  function  the  variational  equation  in  (43)  may 
be  solved  for  the  optical  path,  s,  using  the  Euler-Lagrange  equations. 
At  the  position  of  the  ship  on  the  surface  of  the  ocean  the  incident 
field  will,  appear  to  be  coming  from  the  direction,  f  '(t).  The  error 
angle  between  the  actual  and  apparent  source  directions  is  given  by 
e.  If  we  assume  that  the  refractive  index  is  only  a  function  of  height 
above  the  earth,  with  no  horizontal  gradients,  then  it  can  be  shown, 

[5,  pp.  82-87],  that  the  error  angle,  e,  is  given  by  the  integral 


cot  x  dn 


cot 
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(A3) 


where  ng  is  the  surface  index  of  refraction  and  x  is  the  angle  of  the 
ray  path  with  respect  to  the  horizontal,  as  shown  in  Figure  10.  If  we 
assume  that  n~l  along  the  path,  (since  ngw  1.0003),  the  approximate 
form  on  the  R.H.S.  of  Equation  (44)  is  obtained.  Substituting  the 
refractivity 

N  =  (n  -  1)  x  106  .  (44) 

into  (44)  and  integrating  by  parts  we  obtain 


2 . 4  Time-Domain  Antenna  Response 

In  this  section  we  treat  briefly  the  "voltage"-response*  of  an 
antenna  to  a  source  that  is  moving  slowly  through  its  main  beam  in 
terms  of  a  quasi-static  time-varying  linear  system  approximation.  The 


*When  the  line  into  which  the  antenna  feeds  energy  is  a  waveguide  in 
which  the  dominant  mode  is  one  other  than  a  TEM-mode  the  transmission 
line  voltage  and  current  require  special  definitions.  This  is  well- 
known  and  treated  in  a  number  of  books;  for  example,  see  R.  E.  Collin 
[6,  pp.  145-197]. 


Figure  13  Equivalent  Time-Varying  Linear  System 

/ 

Linear  Transfer 
Operator 
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distant  celestial  source  and  the  moving  antenna  are  shown  in  Figure 

12,  where  the  position  of  .the  source  is  to  be  regarded  as  being  fixed 

in  the  chosen  coordinate  system.  The  unit  vector  p(t)  points  in  the 

direction  of  the  center  (axis)  of  the  main  beam  of  the  antenna  at  a 

given  angular  frequency  of  operation  an  The  unit  vector  p '  is  a  "dummy" 

direction-vector.  The  equivalent  time-varying  linear  system  for  the 

antenna  and  all  sources  of  frequency  in  all  directions  is  shown  dia- 

grammatically  in  Figure  13.  The  driving  term  for  the  time-varying 

linear  system  operator,  L(a>,  t) ,  is  the  Fourier  time-transform  of  all 

incident  electric  fields  of  the  particular  polarization  to  which  the 

antenna  responds;  we  designate  this  quantity  by  E^(p',  uj)  .  The  output 

of  the  system,  V^(m,  t) ,  is  not  a  true  Fourier  transform  because  it 

contains  a  slowly  varying  time-dependence.  The  time-domain  response 

of  the  receiving  system  is  obtained  from  V .  (o>,  t)  by  inverse  Fourier 

A 

transformation. 

vA(t) 


h  /  V“-  c)  ej“t  d“ 


(54) 


The  time-dependent  linear  operator,  L(uj,  t),  involves  an  integration 
over  all  solid  angles  of  the  Fdurier  time-transforms  of  the  incident 
electric  fields  of  particular  polarization  that  the  antenna  is  sensi¬ 
tive  to,  weighted  by  an  "antenna  response  function,"  A(p(t),  p',  ui)  ■ 


VA(o>,  t) 


£p(p', 


ai)  A(p  (t)  ,  p', 


ui)  dft 


(55) 
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This  antenna  response  function  has  a  magnitude  that  Is  proportional  to 
the  square  root  of  the  antenna  power  pattern,  p(p(t),  p',  go)  .  In  the 
case  of  celestial  sources  the  incident  spectrum  ep(p',  go)  will,  of 
course,  be  unknown,  or,  at  best,  will  be  known  statistically  since  the 
processes  that  generate  these  incident  fields  are  random  in  nature. 

If  we  feed  the  antenna  voltage  into  a  narrow  band-pass  R.F. 
amplifier,  as  shown  in  Figure  14,  having  a  transfer  function  as  shown 
in  Figure  15,  output  voltage  of  the  amplifier,  in  the  time-domain,  may 
be  expressed  in  the  form 

vR  R  (t)  =  6 ( t)  [vc(t)  cos  <u0t  -  Vs(t)  sin  Go0tJ  +  Ns(t)  +  NR(t)  (56) 

The  first  term  in  the  R.H.S.  of  (56)  is  the  output  due  to  the  particu¬ 
lar  source  being  observed;  it  consists  of  the  portion  of  the  integral 
In  (55)  over  only  the  extent  of  the  source  being  passed  through  the 
narrow  band-  R.F.  amplifier  with  center  frequency  at  0)o.  The  terms  in 
the  brackets  are  one  of  the  canonical  representation  for  a  narrowband 
signal,  [7,  pp.  132-136],  In  this  case  this  bracketed  term  represents 
the  narrowband  random  signal  due  to  the  source  when  the  antenna  is  not 
moving.  The  3£t)  multiplier  represents  the  very  slowly  changing  ampli¬ 
tude  variation  of  the  received  signal  that  is  due  to  the  motion  of  the 
antenna  pattern.  The  second  term  in  (56),  Ng(t),.is  an  undersired  nar¬ 
rowband  "noise"  term  generated  by  the  extended  integration  in  (55) 
over  all  sources  other  than  that  being  observed.  This  term  should  be 
kept  relatively  small  by  using  a  highly  directive  antenna  with  very 
small  side  lobes  and  by'  observing  a  strong  source  in  a  sparsely  occupied 


Figure  14  Bandpass  R.F.  Amplifier 


Bandpass 
R.F.  Amplifer 


Figure  15  R.F.  Amplifier  Transfer  Function 
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region  of  the  celestial  sphere.  The  tird  terra  in  (56),  N  (t),  is 

thermal  noise  generated  in  the  receiver  which  appears  in  the  output 

of  the  R.F.  amplifier.  This  term  will  depend  upon  the  structure  and 

temperature  of  the  various  stages  of  the  entire  receiver  including 

the  antenna  and  R.F.  amplifier.  The  interferometer  system,  to  be 

discussed  in  the  following  chapter,  incorporates  two  predetection 

antenna/ receiver  subsystems.  The  N  (t)  generated  by  each  of  these 

R 

two  predetection  subsystems  are  statistically  independent  and  hence 
uncorrelated.  The  use  of  the  correlation  detection  system  will  eli¬ 
minate  most  of  the  uncorrelated  thermal  noise  but  will  not  eliminate 

the  correlated,  but  undesired,  noise  terms,  N  (t). 

s 

2.5  Antenna  Tracking  Equations 

In  the  preceding  section  there  appeared  a  multiplicative  term, 

B(t),  that  represented  the  variation  in  the  amplitude  of  the  received 
signal  induced  by  the  changing  position  of  the  source  in  the  power 
pattern  or  "antenna  response  function,"  A(p(t),  p',  to),  of  the  antenna. 
For  a  ship  on  the  ocean  this  relative  motion  of  the  source  in  the  an¬ 
tenna's  power  pattern  is  induced  by  both  the  rotation  of  the  earth 
and  the  motion  of  the  ship  on  the  sea.  If  the  antenna  does  not  track 
on  the  source  to  compensate  for  the  roll,  pitch,  and  yaw  motions  of  the 
ship,  considerable  signal-amplitude  variations  can  result.  These  ampli¬ 
tude  variations  will  be  roughly  periodic  in  nature,  at  least  for  an 
axially  symmetric  main  beam  power  pattern,  and  will  have  a  periodicity 
of  the  order  of  the  mean  value  of  the  roll,  pitch  and  yaw  periods.  The 
magnitude  of  this  amplitude  variation  will  depend  upon  the  ratio  of  the 
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roll,  pitch  and  yaw  angle  variations  to  the  main  beam  width  of  the  an¬ 
tenna. 

In  this  section  we  will  develop  tracking  equations  to  keep  the 
source  in  a  relatively  constant  location  in  the  main  beam  of  the  antenna. 
For  antenna  beamwidths  smaller  than  the  supremum  of  the  standard  devia¬ 
tions,  <5^,  6^,  ^  in  a  given  seastate,  it  will  be  mandatory  for  the 
antennas  of  the  interferometer-navigational-system  to  track  on  the 
source  and  at  least  partially  compensate  for  the  ship's  motion.  The 
equations  which  we  shall  develop  in  this  section  will  be  based  on  as¬ 
sumed  initial  values  of  the  coordinates  of  the  position  of  the  ship. 

This  assumed  position  could  be  based  on  initial  conventional  naviga¬ 
tional  procedures  during  the  startup  of  the  interferometer  or  could 
be  based  on  highly  accurate  position-fixes  obtained  previously  by  the 
interferometer  system. 

The  right  ascension  and  declination  of  a  given  discrete  source 
will  be  known1  in  the  celestial  coordinate  system  of  Figure  7.  The 
cartesian  components  of  this  source  direction  can  be  obtained  in  the 
latitude-longitude  system  of  Figure  6  by  calculating  the  equivalent 
longitude  angle  using  (21)  and  also  using  the  fact  that  6g  =  DEC. 

We  then  use  equations  (16),  (17),  and  (18)  to  obtain  the  cartesian 
components  of  ?(t).  What  we  desire  is  an  expression  for  the  assumed 
(or  actual)  apparent  position  of  the  source  in  the  local  horizon  system 
of  Figure  5  using  an  initial  guess  (or  knowledge)  of  the  position  of 
the  ship  on  the  earth.  If  this  assumed  (or  known)  position  is  given 
by  the  longitude  and  latitude  angles,  A0(t)  and  <5^(0,  respectively, 
we  may  then  obtain  the  assumed  actual  source  position  unit  vector, 
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r^(t),  in  the  local  horizon  system  by  substituting  XQ(t)  and  60(t) 
into  the  A(t)  matrix  in  equation  (24) .  This  gives  us  the  approximate 
(or  actual)  transformation  A0(t),  and  we  obtain 


h'(t)  =  A0(t) . o  r (t) 


(57) 


To  obtain  the  cartesian  components  of  the  assumed  apparent  position 
unit  vector,  ra^(t) ,  we  use  the.  assumed  refraction  matrix,  RQ(t)  , 
given  in  (53),  where,  from  measurements,  the  surface  refractivity  is 
known  and  we  use  assumed  (or  actual)  values  and  A„  .  From  (46) 
we  have 


en(t)  =  10"6N 


+  r 


i  2 

y0 


(58) 


and  from  equations  (49)  and  (50) 


As0(t) 


-  1  1 
tan  1  -y- 


(59) 


The  assumed  apparent  position  unit  vector,  (t) ,  is  then  given  by 


?a0(t)  =  K.0(t)  o  A0(t)  °  r(t) 


(60) 


This  vector  gives  the  direction  in  which  the  optical  axes  of  the  inter¬ 
ferometer  antennas  should  be  pointed.  To  track  the  antennas  in  this 
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direction  we  will  need  equations  of  motion  for  direction  coordinates 

that  are  fixed  with  respect  to  the  ship.  These  coordinates  are  the 

longitudinal,  transverse  and  zenith  coordinates  defined  previously 

for  both  the  "rest"  and  "moving”  ship  positions  as  shown  in  Figure  16. 

The  coordinate  system  x^ ,  x^ ,  xg  is  the  same  as  that  shown  in  Figure  1 

and  moves  with  the  ship.  The  system  x  ,  x  ,  x  defines  the  coor- 

1 0  zo  3o 

dinates  that  the  ship  would  have  on  a  perfectly  calm  sea.  The  roll 
pitch  and  yaw  angles,  f;  (t),  £2(t),  and  E,(t),  are  assumed  to  be  known 
at  all  times  using  the  techniques  described  in  section  2.1. 

What  we  will  now  obtain  are  the  cartesian  components  of  the  as¬ 
sumed  apparent  position  unit  vector  of  the  source  in  the  "moving" 
coordinate  system  (xj ,  x2,  x3) .  By  means  of  this  information  it 
will  be  possible  to  transform  to  any  convenient  angular  coordinate 
system  attached  to  the  ship  (and  the  antenna  mounts),  such  as  a  local 
altitude-azimuth  system,  to  continually  track  on  the  assumed  (or  actual) 
source-position.  By  using  this  scheme  we  shall  get  direct  equations 
of  motion  that  not  only  consider  the  rotation  of  the  earth  but  also  all 
motions  of  the  ship  on  the  sea.  We  obtain  the  final  transformation  in 
two  more  steps.  The  first  of  these  is  a  coordinate  transformation  from 
the  local  horizon  system  to  the  rest-position-cartesian  system.  It  is 
apparent  by  inspection  of  Figure  5  that  this  transformation  is  given 
by  the  matrix 


Tj(t)  = 


cos  A0(t) 


sih  A0(t) 


■sin  A0(t)  0 


cos  A0(t) 


(61) 
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We  now  will  obtain  the  transformation  matrix  from  the  "rest"  system 
to  the  moving  system.  This  can  be  done  in  three  successive  rotations, 
one  for  each  angle  £ .  (t),  5„(t)  and  ?3(t).  The  total  transformation 
matrix  is  the  product  of  these  successive  rotational  transformations. 
Using  Figure  16  the  individual  rotation  transformations  can  be  obtained 
by  inspection.  The  total  transformation  is  thus  found  to  be 


T2(t)  = 


fl  0  0 


0 

COS 

sin 

0 

-sin 

cos 

cos  0  -sin  C. 


0  10 


sin  £2  0  cos  ?2 


cos  sin  £ 


-sin  cos 


0  0 


0 


0 


(62) 
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If  we  denote  the  assumed  apparent  position-unit-vector  as  raQ(t)  in  the 
xx ,  x2 ,  x3  coordinate  system  attached  to  the  ship,  this  vector  can  be 
found  at  any  time  t  from 


rao(t)  =  T2  -  Tt  °  R0  «  Ao  „  f(t)  (63) 

All  calculations  involved  in  this  equation  can  be  done  by  a  com¬ 
puter  and  the  antennas  can  be  servo-controlled  to  track  on  this  posi¬ 
tion.  The  motions  of  thfe  antennas  involved  in  this  tracking  procedure 
will  be  of  a  very  small  angular  velocity  such  as  one-tenth  or  two-tenths 
of  a  degree  per  second.  Such  tracking  motions  should  not  be  difficult 
to  realize.  We  will  assume  that  such  a  tracking  procedure  is  carried 
out  for  each  of  the  two  interferometer  antennas.  Such  a  tracking  pro¬ 
cedure  will,  to  a  very  close  approximation,  make  the  term  B(t)  for  each 
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antenna  a  constant.  This  occurs  because  the  position  of  the  source  in 
the  antenna  power  pattern  remains  fairly  constant  even  though  it  may 
not  be  located  exactly  in  the  middle  of  the  power  pattern  owing  to  an 
initial  slight  error  in  the  assumed  ship  position. 

It  will  be  assumed  later  that  the  term  8(t)  at  each  antenna  of  the 
interferometer  is  never  so  small  as  to  effectively  extinguish  the  source 
signal  when  compared  to  the  background  noise  at  the  output  of  the  corre¬ 
lation  receiver  during  the  course  of  any  measurement  of  position.  The 
time  involved  in  a  single  measurement  is  one  or  two  minutes.  The 
tracking  procedure  must  therefore  "initialize"  the  pointing  of  the  an¬ 
tennas  so  that  the  source  is  well  within  the  main  beams  of  the  antennas 
and  then  keep  it  there  throughout  the  measurement  of  one  or  two  minutes. 
The  initial  pointing,  using  (63) ,  requires  an  a  priori  knowledge  of 
the  ship's  position  on  the  earth  to  within  1/3  or  1/4  of  the  main  beam- 
width  of  the  antennas  in  both  latitude  and  longitude  coordinates.  As 
an  example,  a  dish  3  meters  in  diameter  operating  at  f  =  30  GHz  will 
have  a  main  beam  width,  0j  ~  3.3  x  10-3  radians.  When  observations  of 
the  source  are  made,  close  to  the  zenith,  on  the  basis  of  an  assumed 
ship  position  to  within  about  0,  /3  radians  we  are  guaranteed  that  the 
source  is  initially  located  well  within  the  main  beam  of  the  antenna. 
Thus  \Q  and  6 0  must  be  within  about  1.1  x  10-3  radians  =  .0625  degrees 
of  the  true  ship-position.  This  is  equivalent  to  an  initial  knowledge 
of  position  to  within  about  4.3  miles  of  the  true  position.  Of  course, 
once  the  interferometer  system  is  operating  continuously  the  ship  posi¬ 
tion  will  be  known  at  all  times  to  within  a  distance  of  one-hundredth, 
or  a  smaller  fraction  of  the  above  figure. 
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III.  THE  INTERFEROMETER  NAVIGATIONAL  SYSTEM 

3.1  Time-Varying  Phase  Path  Differences 

In  this  section  we  consider  a  very  general  formulation  for  repre¬ 
senting  the  output  voltages  of  the  two  antenna/R.F.  amplifier  "front 
ends"  of  the  interferometer  in  terms  of  phase  path  differences.  Con¬ 
sider  the  two  antennas  and  R.F.  amplifiers  shown  in  Figure  17.  The 
unit  vectors,  r^(t)  and  s'(t),  are  in  the  direction  of  the  apparent 
source  and  along  the  baseline,  respectively.  The  baseline  distance 
is  given  by  S(t),  and  the  angle  O' (t)  is  the  angle  between  s'(t)  and 
r^(t).  In  vector  notation  we  have 

s'(t)  *  r^Ct)  =  cos  O' (t)  (64) 


As  was  shown  in  section  2.4,  equation  (56),  the  output  voltages  of  the 
first  and  second  bandpass  R.F.  amplifiers  may  be  written  in  the  forms 


VRF<t) 


6 1  ( t)  .[V  <t)  cos  t)0t  -  V3(t)  sin  m0tj  +  Ng(t)  +  NR(t)  (65) 


VRF^t') 


=  62(t)  jyc(t)  cos  (i)Qt  -  vs(t)  sin  w0tj  +  Ng(t)  +  NR(t)  (66) 


The  terms  6^(t)  and  6^(t)  are  slowly  varying  amplitude  functions  gene¬ 
rated  by  the  source  moving  through  the  main  lobes  of  the  antennas.  If 
the  antennas  are  identical  and  are  always  pointing  in  the  same  direction, 
we  have  61(t)  -  82(t).  If  the  antennas  are  tracking  on  the  source  so 
that  the  source  position  does  not  change  in  the  main  beams  of  the  antennas, 


Figure  17  Interferometer  Phase  Path  Geometry 
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^(t)  =  62(t)  =  30»  a  constant  term.  The  terms  in  the  brackets  repre¬ 
sent  narrow-band  random  signals  generated  by  the  incident  fields  from 
the  source.  These  incident  fields  at  each  antenna  differ  by  a  time  de¬ 
lay  factor  due  to  the  different  phase  path  lengths  for  the  Fourier  time 
components  arriving  at  each  antenna.  If  we  assume  that  the  R.F.  ampli¬ 
fiers  have  identical  narrow-band  transfer  functions,  as  shown  in  Figure 
15,  and  that  the  frequency  dependence  of  the  index  of  refraction  profile, 
n(£  ',  w) ,  is  smooth  in  the  frequency  neighborhood  of  u0,  then,  to  a 
close  approximation,  each  Fourier  component  of  the  incident  random, 
narrow-band  source-signal  will  have  the  same  time  delay.  We  may  then 
write 

Vc(t)  cos  oaQt  -  Vs(t)  sin  w0t  =  Vc(t  -  x)  cos  m0(t 

sin  u>0 ( t  -  t) 

The  time  "delay,"  i(t),  is  a  function  of  time  given  by  the  difference 
between  the  phase  path  integrals  from  the  source  to  the  second  and 
first  antennas,  respectively,  divided  by  the  free  space  velocity  of 
light,  c.  Note  that  a  negative  "delay"  corresponds  to  a  time  advance¬ 
ment  of  the  signal  at  antenna  if 2  compared  to  the  signal  at  antenna  ffl. 


-  i)  -  V  ( t  -  x) 
■  S1 

(67) 


tu>  -  t 


* 

f  n(T^  (t)  ,  m0)  d&'2  -  J  nU'(t),  oj0)  dJlJ 


(68) 


P£t) 


Px(t) 


Where  P^,  and  are  the  actual  phase  paths  from  the  source  to  the  second 
and  first  antenna,  as  shown  in  Figure  17,  and  0>o)-  and  n(£^(t),  to  ) 
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are  the  respective  indices  of  refraction  at  the  center  frequency,  w  , 
along  the  phase  paths.  We  now  replace  the  indices  of  refraction  in 
equation  (68)  by  their  representations  in  terms  of  ref ractivity,  N, 
using 


n  =  1  +  N  x  10~6 


(69) 


and  note  that  to  a  very  good  approximation 


/ 


d 4-  d(t) 


(70) 


where  d(t)  is  the  straight  line  approximation  to  the  last  section  of 
phase  path,  ?2,  as  shown  in  Figure  17.  The  result  is 


x(t) 


Looking  at  Figure  17  we  see  that 


d(t)  =  S (t)  cos  ©' (t) 


S(t)  §'(t)  •  ?;(t) 


(72) 


which  yields  the  formula 


T(t) 


S(t) 


c 


§' (t)  *  fg  (t)  +  Ax  (t) 


(73) 
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where  Al(t)  is  the  second  term  on  the  R.H.S.  of  (71).  If  we  assume  that 
of  the  refractivity  of  the  atmosphere  is  only  a  function  of  height  (ra¬ 
dially  stratified  model),  the  integrals  in  Ar(t)  over  P^  and  -  d(t) 
are  equal.  If,  then-,  we  approximate  N(l^,  w0),  in  the  last  section  of 
integration  in  P^,  of  length  d(t),  by  Ng(t,  to0),  the  refractivity  at  the 
surface  which  is  also,  in  general,  a  function  of  time,  we  obtain 

At  ( t)  =  Ns(tj  Uq)  d(t)  (74) 

Substitution  into  (73)  and  the  use  of  (69)  and  (72)  yields 

ns(t,  u)  )  S  (t) 

r(t)  *  - - -  s'(t)  •  f'(t)  (75) 

C  a 

Let  us  define  the  time-varying  phase  difference 

ns(t,  oj0)  S  (t)a;0 

cjj(t)  =  uJqT  ( t )  =  — — — - 1 -  s'(t)  •  r'a(t)  (76) 

c 

We  may  then  rewrite  (67)  as  ■ 

V  (t)  cos  w0t  -  V  (t)  sin  m0t  = 

2  s2 

VC](t  -  t)  cos  (m0t  -  <j>(t))  -  Vg^t  -  t)  sin  (u0 1  -  ^ ( t ))  (77) 

The  noise  terms  Ne  (t)  and'  Nc  (t)  are  generated,  as  explained  in 
S1  b  2 

section  2.4,  by  sources  other  than  that  observed.  These  two  noise 
terms  will  be  weakly  correlated  but  will  not  differ  by  a  simple  time 
delay  since  they  are  each  generated  by  fields  incident  from  sources 
at  all  angles  weighted  by  the  antenna  response  function,  as  shown  in 
equation  (55) . 
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The  receiver  noise  terms,  NR  (t)  and  (t) ,  are,  in  general,  un- 

1  2 

correlated  and  their  effects  are  greatly  reduced  by  the  correlation  re¬ 
ceiver  to  be  described  in  the  following  section. 

3.2  The  Correlation  Receiver 

We  proceed  now  to  consider  the  operation  of  a  basic  correlation 

i 

receiver  to  be  used  for  navigation  by  observing  a  discrete  radio  source 
whose  celestial  coordinates  are  known  as  a  function  of  time.  In  the 
following  section  we  will  present  the  technique  for  utilizing  the  in¬ 
terferometer  and  correlation  receiver  for  navigation.  The  basic  cor¬ 
relation  receiver  is  shown  in  block  diagram  form  in  Figure  18.  The 
output  signals  from  the  bandpass  R.F.  amplifiers  are  of  the  form  given 
in  the  preceding  section  by  equations  (65)  and  (66) . 

vRF  (t)  =  62(t)  jVc(t)  cos  t)0t 

VRF  ^  =  pc^)  cos  wot 

2  L  2 

For  the  interferometer  navigational  system  to  be  described,  we  as¬ 
sume  that  the  identical  antennas  making  up  the  elements  of  the  interfero¬ 
meter  are  tracking  on  the  assumed  source  position  so  closely  that  the 
source  stays  well  within  the  main  beams.  Therefore,  the  terras  6^(t)  = 
62(t)  =  B0(t)  change  very,  very  slowly  and  never  become  so  small  that 
the  signal  from  the  source  gets  lost  in  the  noise  at  the  output  of  the 

correlation  receiver.  We  assume  also  that  the  R.F.  and  I.F.  amplifiers 

i 

of  both  predetection  subsystems  have  respectively  identical  bandpass 


-  Vs(t)  sin  m0tj  +  Ns(t)  +  NR(t)  (65) 

-  V„(t)  sin  a)  tl  +  Ns(t)  +  NR(t)  (66) 

2  °  J  2  z 


Figure  18  The  Interferometer  and  Correlation  Receiver 
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characteristics  such  as  are  shown  in  Figure  19  (p,  68).  The  R.F.  am¬ 
plifiers  are  "narrow-band";  the  refractivity  function,  N ( 2. '  ,  w) ,  u> 
assumed  to  be  smooth  with  respect  to  frequency  in  the  neighborhood  of 

w  .  We  can  thus  represent  the  relative  variations  of  the  two  source- 
o 

signals  v  and  v  by  a  simple  time  delay  (or  advancement)  as 
Hr  ^  Rt  g  1 

explained  in  the  previous  section.  Using  equations  (65),  (66)  and 
(77) ,  we  have 

w0t  -  Vg(t)  sin  u>0 1 J  +  Ns(t)  +  NR(t)  (78) 

cos  poQt  -  4>(t))  -  Vs(t  -  t)  sin  (w0t  -  4>  ( t 

(79) 

where  from  (75)  and  (76) 

ng(t,  mQ)  S (t) 

t (t)  =  - - -  s'(t)  •  S^(t)  (75) 

w  ns(t,  m0)  S (t) 

$(t)  = - -i -  s’(t)  •  f’(t)  (76) 

c  a 

The  signals  (78)  and  (79)  are  mixed  with  signals  from  the  local  oscil¬ 
lators  as  shown  in  Figure  18.  These  local  oscillators  are  phase  con¬ 
trolled  with  voltages  given  by 


v,10 


:t)  p 


=  B0(t)  fVr(t)  COS 


VRF  (t)  -  6o(t)  [vc(t  -  T) 


+  Ns(c)  +  %(t) 
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vLo  (t)  =  2  cos  jjCoJo  "  'J>IF)  tj  (80) 

vLo  (t)  =  2  cos  J(aj0  -  uIF)  t  -  4>0(t)J  (81) 

The  phase  term  <j)0(t)  will  be  explained  in  the  following  section.  The 
mixed  signals  are  given  by 

v  (t)  =  6 0 C t )  (Vc(t).  cos  a)  t  -  V  (t)  sin  cjt  t  +  f"Ns(t)  +  NR(t)lv  (t) 

mi  Li  rt  °i  k  j  \  L°i 

+  fl  ((2u0  -  ^IF)  t)  (82) 


and  and  f  are  the  upper  side  band  terms  at  2oo0  -  that  will  be 

eliminated  when  the  signals  pass  through  the  I.F.  filters.  The  signals 
at  the  I.F.  outputs  will  have  the  forms 
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vip(t)  =  G0B0  Vc(t)  cos  U];Ft  -  Vg(t)  sin  +  N^Ct)  (85) 

vIf,(t)  =  GO0O  Vc(t  -  t)  cos  (u)IFt  +  A<Ji(t))  -  Vs(t  -  x)  sin  (o)Ipt  +  A C t ))J 

+  NTTJ(t)  (86) 

IFz 

where  G^  is  the  I.F.  amplifier  gain  and  NjF(t)  and  Njp(t)  are  the  total 
non-signal  ncise  terms  at  the  outputs  of  respective  I.F.  amplifiers. 

The  time  delays  at  the  outputs  of  the  I.F.  amplifiers  are  con¬ 
trolled  together  to  recorrelate  the  terms  Vc(t  -  x)  and  Vg(t  -  x)  in 
(86)  with  their  non-delayed  counterparts  in  (85) .  This  can  be  done  in 
one  of  two  ways  depending  upon  the  sign  of  the  "delay,"  x(t).  The 
sign  of  this  term  will  be  known  even  though  its  exact  magnitude  will 
be  unknown  when  the  ship's  exact  position  is  unknown.  For  t  >  0  we  set 
x„  =  0  and  set  ij(t)  to  an  assumed  value  for  x(t)  based  on  the  assumed 
ship-position.  This  is  done  using  the  assumed  unit  source  direction 
vector,  raQ(t),  given  in  equation  (63)  of  section  2.5.  Denoting  the 
assumed  apparent  position  of  the  source  in  the  local  horizon  frame  by 
the  unit  vector.,  ra'(t),  we  define,  using  (75), 

A  ^ S(t) 

T0(t)  =  — — - — —  §’(o  •  ta;( o  (87) 

Now,  using  the  invariance  of  the  scalar  product  and  the  fact  that  in 
the  "moving"  system  attached  to  the  ship  the  unit  baseline  vector  has 
the  representation  X. ,  we  obtain  from  (87) 
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T0(t)  = 


ns(t,  u)D)  S  (t) 


(88) 


where  the  X  component  of  raQ(t)  will  be  known  using  (64).  Since 

V  (t)  and  V  (t)  are  random  signals  of  bandwidth  Av,-,  their  auto 
i  53  ]_  lr 

correlation  coefficients  approach  zero  very  rapidly  for  time  separations 
greater  than  — ^ —  .  It  is,  therefore,  mandatory  that  the  error  in 


IF 


T0(t)  be  only  a  fraction  of 


Av 


Substituting  ns  w  1  and  denoting 


IF 


the  angles  from  the  s’  vector  to  the  actual  and  assumed  apparent  unit 
direction  vectors  r^  and  as  0^  and  0a^,  respectively,  we  obtain  with 
[rao]x  =  cos  e;o 


AX 


At  =  t ( t)  -  Tc(t)  «  — 


cos  Ol  -  cos  Gl  << 


IF 


IF 


S(t) 


(89) 


where  AX  is  the  equivalent  wavelength  of  the  I.F.  bandwidth  fre- 
1  r 

quency.  For  observations  near  the  zenith  we  expand  the  cosine  func¬ 
tions  in  Taylor's  series  about  0  =  t/2.  The  result  is 


AX 


a©  =  el  -  e,,  << 

^o 


IF 


a„  -a  "  S(t) 


(zenith  observations)  (90) 


If  we  assume.  Av  =  200  MHz  and  S(t)  =  200  meters,  we  find,  for  example, 
-L  F 


A0  <<  7,5  x  10  3  radians  =  26  min 


(91) 


Requirements  for  accuracy  of  tb)is  order  of  magnitude  will  not  introduce 
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additional  problems  in.  the  procedure  since  use  of  the  antenna  tracking 
technique  in  section  2.5  requires  a  A@max>  that  is,  a  fraction  of  the 
main  beam  widch.  Such  antenna  tracking  requirements  will  usually  satisfy 
the  additional  requirement  in  (89)  easily. 

For  the  case  of  x  <  0  we  set  x^  =  0  and  set  x2  =  T0(t) .  The  time 
delay  circuits  in  Figure  18  may  consist  of  binary  sequenced  delay  lines 
such  as  are  now  used  in  the  Hat  Creek  interferometer  system.  In  using 
such  a  discrete  delay  circuit  the  resolution  must  be  kept  small  enough 

so  that  the  error,  Ax,  is  always  less  than  1/A\>  c. 

IF 

The  inputs  into  the  final  multiplier,  v^t)  and  v(t),  as  shown  in 
Figure  18,  will  have  alternate  forms,  depending  upon  whether  x(t)  is 
positive  or  negative.  These  forms  are  given,  for  x  >  0  and  x  <  0, 
respectively,  by  x  >  0  =t>  x0  =  0,  x  =  xQ  =  x  -  Ax 


Vj  (t) 


G06(t 


-  xG)  Vc(t  -  t  +  Ax)  cos  w1F(t  -  tq) 


-  Vs(t  -  x  +  Ax)  sin  oiIF(t  -  xQ) 


J  +  NIF(t  -  x 0 


)  (92) 


v2(t) 


G0S(t) 


Vr(t  -  x)  cos  (mT_t,  -  A<J)> 

L  1  ■■  IF 


-  Vs(t  -  t)  sin  (uijpt  -  Ad)J  +  Nip (t) 


(9  3) 


x  <  0  =>  x  =  G,  x2  =  x0  =  x  -  Ax 


vT  (t) 


G03(t) 


cp'  UJip,_  's£-/  “'ip1-  +  NIF^ 


V r ( t )  cos  uTT?t  -  Ve(t)  sin  wT1?t 


(t)  (94) 
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(t)  =  G0fJ(t  -  TQ)jvc(t  +  At)  cos  (wIp(t  -  Tq)  - 


-  Vg(t  +  At)  sin  (ii>IF(t  -  tq)  -  A^j  +  NjF( 


t  -  t0) 


(95) 


The  paired  inputs  in  (92)  and  (93)  or  (94)  and  (95)  are  now  mul¬ 
tiplied  and  the  resultant  put  through  the  final  low-pass  filter  or 
"integrator."  The  transfer  function  of  this  low-pass  filter  is  shown 
in  Figure  19.  The  products  put  into  this  filter  are  of  the  forms 

T  >  0  ; 


Vj(t)v2(t)  =  G£8(t)3(t  -  t0) 


Vc(t  -  i  +  At)  vL(t  -  t)  +  V_(t  -  t  +  At)  Vs(t  -  t) 
1  1  1  1 


cos  ( A 4>  -  lo^Tq) 


+ 


Vc(t  -  t  +  At)  Vs(t  -  t)  -  Vc( t 


-  t)  Vo(t  —  T  +  At) 


sin  (A<f>  -  miFT0)} 


+  g(2uipt)  +  N0(t) 


(96) 


<  0 


Vi(t)v2(t)  =  (t)  8  (t  -  tc)  >: 


Vc(t  +  At)  V c ( t )  +  Vs(t  +  At) Vs  ( t) 
1  1  1 


cos  (A cp  +  dqift0!)) 


+ 


V_  (t  +  At  )  Vs  (t)  -  Vc(t)Vg(t  +  AT) 
C1  1  ■  1  1 


sin^Ad  +  Cj0xFTO)|  +  h  (Zui^t)  +  N0(t)  (97) 
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The  terms  g  and  h  represent  upper  side  band  terms  at  a  frequency  of 
that  will  be  eliminated  by  the  low-pass  filter.  The  terms  N0 

and  Nn  are  noise  terms  with  bandwidths  of  2m,.  _.  When  the  product 

2  IF  -  L 

in  (96)  or  (97)  is  passed  through  the  low  pass  filter  only  the  low 
frequency  components  of  the  noise  terms  will  remain.  The  time  average 
noise  power  is  directly  proportional  to  the  noise  bandwidth  so  we  must 
make  the  low  pass  cutoff  frequency  W  as  small  as  possible  while 

Lr 

still  passing  the  first  term  in  (96)  or  (97).  Note  also  that  since 
At  is  much  smaller  than  the  reciprocal  of  the  bandwidth  of  Vc  or  Vs 
we  will  have  with  t'  =  t  or  t'  =  t  -  x 


Vc(t'  +  Ax)V„(t')  +  Vs(t'  +  Ax)Vo(t') 
.1  1  1  L  1 


v§(t’) 


vl(t') 


(98) 


and  . 


Vc(t'  +  Ax) Vs ( t ' ) 
i  1 


Vc(t’)VB(t'  +  Ax) 
C1  S1 


0 


(99) 


The  final  output,  will  have  the  form 

) 

vLp(t)  =  K(t)  [<V*(t')>  +<v|(t')>]  cos  (A<p  ±  wIFx0)  +  n0(t)  (100) 


where  K  >  indicates  the  passage  of  the  enclosed  term  through  the  low 
pass  filter.  If  this  low  pass  filter  were  an  ideal  integrator  these 
brackets  would  indicate  time  averaging.  The  term  K(t)  is  a  very  slowly 
varying,  almost  constant,  function  whose  origin  lies  in  the  antenna 
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amplitude  function  B(t).  The  term  n0(t)  is  a  very  small  noise  term 
residual  that  was  passed  through  the  low  pass  filter.  Typical  sample 

functions  for  the  terms  K(t),  [  v£(t')  +  v|(t')  ],  and  nn(t)  are  shown 

1  S1  1 

in  Figure  20.  For  the  sake  of  brevity  we  will  denote  the  multiplicative 
term  of  the  cosine  function  as  V0(t). 


V0(t)  =  K(t)[<v2(tr)>  -j-  < V§  ( t '  )> 


(101) 


This  term  will,  in  general,  be  almost  constant  with  a  slight  quasi¬ 
harmonic  variation  about  the  mean  value  <VQ>  as  shown  in  Figure  21. 

The  minimum  period  of  this  quasi-harmonic  variation  will  be  limited 
by  the  low  pass  filter  and  will  be  approximately  equal  to  the  reci¬ 
procal  of  v  .  The  minimum  value  of  W  will  be  set  by  the  requirement 

Lji  Li 

that  we  must  pass  the  cosine  term  shown  in  (96)  or  (97).  The  argument 
of  this  cosine  term  is  the  key  .to  using  the  interferometer  for  navi¬ 
gation.  To  pass  this  cosine  term  the  low  pass  filter  must  have  an 
upper  cutoff  frequency  of  at  least 


LP 


mm 


d_ 

dt 


(^4=  ( t)  ±  C0IFTo(t)) 


max 


(102) 


We  will  now  define. the,  local  oscillator  phase  shift  term  <t)0(t)  an^ 

1 

show  how  the  interferometer  may  be  used  for  navigation.  For  a  more 
comprehensive  discussion  of  the  correlation-receiver  and  other  types 
of  receivers  the  reader  is  referred  to  such  references  as.  Kraus, 

[3,  pp.  236  -  290] ,  Bracewell,  [8],  or  Christiansen  and  Hogbom,  [9, 
pp.  190  -  210]. 


Relative  Amplitude 
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Figure  21  The  Amplitude  Multiplier  of  the  Cosine  Function 


53 


3.3  The  Navigational  Technique 

All  work  up  to  this  point  was  done  to  lay  a  foundation  lor  this 
section.  As  was  shown  in  the  preceding  section  the  output  of  the 
correlation  receiver  with  the  interferometer  antennas  tracking  on  a 
discrete  source  is  given  by 

vLp(t)  =  V0(t)  cos  A 4> ( t )  ±  iDIFT0(t)  +  nQ(t)  (103) 

'  / 

The  phase  difference  term  is  given  by 

A<{>(t)  =  -  4>0(t)  ,  (84) 

where  the  term  <J>(t)  is  the  unknown  fringe  phase  of  the  source  and 

<f>0(t)  is  the  known  phase  injected  into  the  local  oscillator  it 2  as 

shown  in  Figure  18.  The  method  to  be  presented  in  this  section  will, 

in  principle,  work  for  almost  any  4>0(t)  that  is  chosen.  There  is, 

however,  an  optimum  choice  for  cf>0(t).  Recall  that  in  the  last  part 

of  the  previous  section  we  discussed  the  minimum  low  pass  filter  cutoff 

frequency.  Since  the  time-average  noise  power  in  n  (t)  is  directly 

proportional  to  W  we  should  make  this  cutoff  frequency  as  small  as 
Lr  , 

possible  while  still  passing  the  cosine  term.  The  minimum  cutoff  fre¬ 
quency  required  to  accomplish,  the  passage  of  the  cosine  term  is  given 
in  equation  (102) .  It  is  apparent  that  a  minimum  W  will  be  allowed 
by  minimizing  ( t )  at  all  times.  The  optimum  choice  for  bo^t)  is  thus 
the  best  guess  we  can  come  up  with  for  <j)(t).  This  best  guess  is  obtained 


54 


directly  from  the  antenna  tracking  procedures  developed  in  section  2.5, 
as  was  also  the  best  guess  for  tQ(t)  that  was  obtained  in  the  previous 
section.  If  we  obtain  the  best  guess  of  ra  (t)  in  the  "moving"  coor- 
dinate  system,  (x^,  x2 ,  x^),  attached  to  the  ship  using  equation  (63) 

t  aQ ( t )  '  T2(t)  •  T  j (t)  •  R0(t)  •  A0(t)  •  r ( t) ,  (63) 


where  T^,  T^,  R0  and  A0  are  given  using  the  procedures  in  section  2.5, 
the  best  guess  for  $(t),  namely  <j>Q(t)  ,  will  be  given  by 


4o(t) 


wcns ( t) S ( t) 
c' 


(104) 


This  is,  of  course,  just  equal  to  w0x0(t).  The  rationale  behind  the 
"best  choice"  for  <}>0(t)  can  be  summarized  by  stating  that  we  wish  to 
make  the  "frequency  aperture"  of  the  low  pass  filter  as  small  as  pos¬ 
sible  to  restrict  the  passage  of  the  noise  term.  At  the  same  time  we 
need  to  collapse  the  frequency  spectrum  of  the  phase  term  <£(t)  so  it 
will  fit  through  the  low  pass  frequency  aperture.  To  accomplish  this 
we  subtract  a  best  guess  for  <|)(t),  namely  4>0(t).  Once  the  collapsed 
phase  term  passes  through  the  low  frequency  aperture  we  add  back  the 
known  term  $0(t)  and  obtain  <f>0(t)  with  a  minimum  noise  transfer. 

We  wTill  now  discuss  how7  we  add  back  the  term  $G(t)  onto  the  output 

phase.  We  will  also  subtract  or  add  the  known  term  w  T0(t)  so  we  obtain 

lr 

the  term  <Kt)  alone.  The  correlation  receiver  output,  v^p(t)  is  sketched  in 
Figure  22  including  the  noise  term  nQ(t)  and  the  amplitude  term  VD(t). 

Since  the  low  pass  filter  is  chosen  to  just  pass  the  cosine  term  the  periods 
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of  variation  of  V0(t),  n0(t)  and  cos  (A<j>  ±  o)^pT0)  will  be  about  the 

same.  The  noise  term  nQ(t)  and  the  amplitude  term  VQ(t)  will  inherently 

introduce  slight  errors  in  the  measurement  of  the  phase  of  the  cosine 

term.  The  output,  can  be  into  a  phase  detector  directly 

or  can  be  prepared  for  phse  detection  using  a  shaping  circuit.  The 

phase  detector  can  not  detect  the  absolute  phase  of  the  cosine  term  but 

will  give  the  relative  changes  of  phase  of  the  cosine  term.  There  is 

thus  an  ambiguity  of  2Ntt  in  the  measured  phase  of  v  (t).  This  does 

Li 

not  matter  since  we  will  only  be  interested  in  phase  variations  of 
<J>(t).  To  obtain  these  phase  variations  of  c(j ( t )  we  use  the  system  shown 
in  Figure  23.  The  input  to  the  phase  detector  is  the  correlation-re¬ 
ceiver's  output,  v^p(t),  in  analog  form.  The  phase  detector  gives  a 
digital  output  of  the  term 


$(tK)  =  A<j)(tK)  ±  iuIFT0(tK)  +  2Ntt . 


(105) 


where  2Ntt  is  the  unknown  ambiguity  in  the  absolute  value  and  t^  are 
the  discrete  time  values  for  which  this  function  is  read  out.  The 
resolution  of  tj^  must  be  small  enough  to  track  the  variation  of  3>(t) 
closely.  The  measured  numerical  phase  is  then  read  into  the  adder 
circuit  which  adds  the  known  <j)0(tK)  and  ±  wIFT0(tjr)  to  $(tK).  The 
numerical  output  of  the  adder  will  be  +  2Ntt.  Using  this  output 

we  can  obtain  the  change  of  phase  of  $(t)  between  any  two  discrete  times 
t^  and  tj  by  simple  subtraction.  Let  us  define  this  measured  phase 


difference  as 
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Figure  22  Correlation  Receiver  Output  (Dashed  lined  indicates 
the  output  without  noise  and  with  a  constant 
multiplier  Ve  ) 


Vl.p.  (t ) 


Figure  23  Fringe  Phase  Recovery  System 


vL.p.(t)  — H 


Adder 


+aJIFrO^V 
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Y  <t±  |  tj)  =  4>  ( t±)  -  *<tj)  (106) 

To  see  how  we  will  use  this  measured  phase  variation  to  obtain 
the  position  of  the  ship  we  use  (57)  and  (76)  to  write  d> ( t ±)  as 

w0  =  = 

Kti)  =  ~  ns(tit  to0)  s(t±)  ^'(ti)  •  R(ti)  •  A(t±)  •  r(ti)  (107) 

Using  (106)  we  obtain 

wo  = 

Y (ti  I  tj)  =  —  ns(t1)S(t1)s' (ti)  •  R(t±)  ■  A(ti)  *  r(t±) 

~ns (t j ) S (t j ) s ' (t j )  •  R(tj)  *  A(tj)  •  ?(tj)  (108) 

A  few  approximations  can  be  made  for  this  expression.  The  first  of 
these  is  the  assumption  that  ns(t)  is  approximately  constant  over  the 
course  of  several  minutes.  The  second  approximation  will  be  to  replace 
R  with  the  known  matrix  RQ.  These  approximations  are  made  only  to  sim¬ 
plify  the  navigational  technique  and  should  introduce  no  serious  errors. 
The  success  of  the  method  is  not  contingent  upon  these  assumptions  and 
a  more  general,  but  messy,  technique  can  easily  be  developed  along  the 
same  lines  as  our  simplified  method,  but  without  the  above  assumptions. 
Substituting  the  above  assumptions  into  (108)  we  obtain 

y(ti  tj)  =  — ■  -S  ,  S(t1)§'(ti)  •  R0(ti)  *  A(ti)  •  £(ti) 

-  S(tj)8'(t  )  •  R0(tj)  "  A(tj)  *  b(tj ) 


(109) 
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When  we  look  at  this  expression  we  notice  that  .all  terms  on  both  the 
L.H.S.  and  R.H.S.  of  this  equation  are  known  or  measured  except  A(t^) 
and  A  (t^).  That  is 

a)  y ( t^ | t j )  is  measured. 

b)  S(t),  the  baseline  length,  is  measured. 

c)  s(t),  the  unit  baseline  vector  is  measured  in  the  hori¬ 
zon  system.  Its  cartesian  components  are  given  in  equa¬ 
tions  (11)  to  (13). 

d)  Rq ( t)  is  known  using  the  technique  of  section  2.4. 

e)  r(t),  the  unit  source  direction  vector  in  latitude-longi¬ 
tude  coordinates,  is  known  from  tabulated  data. 


To  simplify  the  notation  to  follow  let  us  define  the  vector 

_  w  n.  _ 

T(t)  =  S(t)  s’(t)  R  (t)  (110) 

c  o 

The  components  of  this  vector  will  be  known.  Substituting  into  (109)  we 
obtain  the  equation 

B(ti)  «  A(t±)  •  r(ti)  -  B(tj)  «  A(t^)  r(tj)  =  Y (t^l  t ^ )  (111) 

Using  the  1,2,3  representation  for  the  cartesian  components  of  B  and 
r  we  can  rewrite  equation  (111)  in  the  form 


££  VV  W  'l11!1  •  VV  amJ(tj)r»(tj)  ‘1'(ti|tj> 

m=!  £=1  (112) 


On  examining  A(t)  in  equation  (24)  we  notice  that  its  unknown  elements 
a^^(t),  if  found,  would  give  us  the  position-longitude  and  latitude  co¬ 
ordinates  of  the  ship  on  the  earth.  Let  us  say  we  desire  to  know  the 
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position  of  the  ship  at  t^  =  tQ.  A  knowledge  of  the  matrix  elements 
a^(to)  will  provide  us  with  the  ship's  position.  We  notice  in  equa¬ 
tion  (112)  that  there  are  two  sets  of  unknowns,  a.  (t.)  =  a.  (t  )  and 

xm  j  xm  o 

a„  (t. ).  To  express  a.  (t.)  in  terms  of  the  desired  unknowns  a,  (t  ), 

Xm  i  xm  i  xm  o 

where  t^  >  t^,  we  use  the  equations  of  evolution  derived  in  section  2.2. 

Using  these  equations,  (34)  to  (41),  we  substitute  for  a  (t^)  in  terms 

of  a.  (t  )  in  (112).  The  result  can  be  expressed  in  the  simplified 
Xm  o 

form 

3 

E 

m=l  X— 1 

Where  the  coefficients  F  are  given  on  the  following  page  in  terms 

uLv 

of  known  or  measured  functions.  The  functions  f^,  g,  and  are 
given  in  equations  (29)  to  (32)  and  are  measurable.  Using  (25),  the 
matrix  A(tQ)  is  given  by, 

_  cos  A(t  )  -sin  \(t  )  0 

o  o 

A(tQ)  =  sin  X(tQ)  sin  <5 ( t cos  X(tQ)sin  6(tQ)  -cos5(tQ) 

-  sin  \(t  )  cos  6(t  )  cos  \(t  )cos  6(t  )  sinS(t  ) 

o  o  o  o  o 

(114) 

To  obtain  f^(t)  and  f^(t)  from  (29)  and  (30)  we  use  the  assumed  value 
60 (t)  for  the  declination  6(t). 


(ti'  ^ 


a_,  (t  )  =  y  (t .  t  ) 

mil  o  i  o 


(113) 
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The  values  of  F  .  are: 

td£ 


11 


12 


-  <co> 


(t±  lto)  =  Bi(ti)j^f1(ti)r1(ti)  -  f2(t.)r2(t.) 

(tilt:o)  =  B1(ti)^f2(t1)r1(t1)  +  f1(ti)r2(t.)J  (tQ) 


F21^Ci  ,to)  ~  B2^ti') 


+  B. 


-  B2(tQ)r1(t0) 


’22 


=  B2(ti)  [^2  8i  (t±)  Ti  (ti>  + 

+  B3(tI)  ^2(ti)g2(t±)ri(ti)  +  f1(ti)62(ti)r2(ti)] 


-  B2(to)r2<to) 


F23(ti!to)  =  B2(tl)  [8l<ti)r3(ti)  +  g2(ti)r3(ti)]  •  B2(tO)r3(tO) 
=  B2(ti)  [fl(ti)g2(ti)ri(ti)  +  f2(ti)82(ti)r2(ti)] 
,(ti} 


+  B 


-  B  (t  )r  (t  ) 
3  o  1  o 


F32(ti'tO)  =  ~W 


+  B  (t.) 

3  1 


f2(ti)g2(ti)ri(ti)  +  f1(ti)g2(ti)r2(t.)J 

_  ^2 ( fci) ® l <tl>rl<ti>  +  £l<ti)S1(tl)l2(ti)] 


-  B  (t  ) r  (t  ) 
3  o  2 '  o 
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"  r3^ti)  f B3^Ci)8i^ti^  ~  B2^tI)82(ti)|  "  B3(tO)r3(tO) 


(115) 
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There  are  several  algebraic  techniques  for  extracting  the  position 

coordinates  X(c  )  and  6 ( t  )  from  the  set  of  equations  in  (113).  All 
o  o 

of  these  techniques,  but  one,  require  the  solution  of  a  nonlinear  and 
mixed  system  of  equations.  Looking  at  (114)  we  see  that  there  are 
only  two  "independent"  variables,  A(tQ)  and  8(tQ),  that  generate  the 
eight  nonzero  elements  of  the  A  matrix.  These  independent  variables 
generate  the  matrix  elements  in  a  nonlinear  manner  through  transcendental 
equations.  Even  though  there  are  only  two  "independent"  variables  X(tQ) 
and  S(tQ)  the  equation  (114)  can  generate  eight  "linearly  independent" 
equations  in  the  eight  linearly  independent  variables  am£(t0)  (with 
a13=  0  ).  As  an  example  of  the  linearization  of  a  nonlinear  system  con¬ 
sider  the  equation 


sin  XQ  +  y  (t)  cos  xQ  =  Z(t) 


(116) 


This  is  a  nonlinear  transcendental  equation  in  the  single  unknown  XQ* 
If.  we  let  Wq  =  sin  xq  then  this  becomes  a  nonlinear  equation  in 


W 


o 


+  y(t) 


-  Z(t) 


(117) 


The  functions  y(t)  and  Z(t)  are  parametric  in  the  variable  t.  They 
are  not  independent  but  can  assume  different  but  coupled  values  as  t 
is  changed.  The  system  (one  equation)  in  (117)  can  be  solved  for  Wq 
directly  using  only  one  value  of  t.  We  can,  however,  also  solve  the 
system  by  defining  a  second  linearly  independent  variable  =  cosxQ 
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and  rewriting  (116)  as 


Wo  +  y(t)  Wx  =  Z(t) 


(118) 


If  we  let  t  take  on  the  values  t^  and  we  obtain  the  coupled  system 
of  two  linearly  independent  equations 


W 

o 


+  y(t  )  w  = 

a  ■  J. 


z  (t  ) 

a 


wo +  y(V  wi =  z  (tb} 


(119) 


which  can  be  solved  directly  by  elimination. 

In  the  same  manner  as  we  converted  the  nonlinear  equation  in  (116) 

to  the  linear  system  in  (119)  we  will  now  convert  the  nonlinear  system 

in  the  two  variables,  A(tQ)  and  6(tQ),  in  (113)  to  a  higher  order  system 

in  the  eight  linearly  independent  variables  a^(to).  Thus,  the 

technique  for  obtaining  the  ship's  position  at  any  time,  t^,  is  to 

measure  and  calculate  the  values  F  . (t .  t  )  and  y(t.  t  )  for  eight 

m£  i  o  i  o 

consecutive  discrete  times  t^  >  tQ.  We  then  simply  solve  the  linear 
system  of  eight  equations  in  eight  unknowns. 


3  3 


F  (t  It  )  a  (t  )  =  y(t  It  )  i  =  1,1 

m=]_  £=1  mx  l  o  mx  o  i  o 


(120) 


We  need  only  to  solve  this  linear  system  for  two  unknowns  such  as 


and 


a  (t  )  =  cosA(t  ) 
11  o  o 


a  (t  )  =  sin5 (t  ) 
3a  o  o 


(121) 


(122) 
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to  obtain  the  ship’s  position  at  any  time,  t  ,  This  solution  can 
be  obtained  easily  using  numerical  methods  such  as  Gauss'  elimination 
or  even  a  direct  use  of  Cramer's  rule.  Note  that  our  solution  for 
position  will  be  at  a  previous  time,  tQ,  To  obtain  the  position 
measurement  at  a  future  time  we  obtain 

a  (t  )  =  -Jl  -  a  2  (t  ) 

1 2  0  T  11  o 

and 

a  (t  )  ■  -"vf  1  -  a  2  (t  ) 

23  o  *  33  o' 

then  directly  use  the  equations  of  evolution  in  (35)  and  (42) 

cos  A  (t)  =  f  (t)  a  (t  )  +  f  (t)  a  (t  )  (125) 

1  11®  2  12® 

sin  <5  ( t )  =  g  (t)  a  (t  )  +  g  (t)  a  (t  )  (126) 

1  33  °  2  23  ° 

In  the  preceding  sections  we  used  the  term  "assumed  (or  known)" 
ship  position,  unit  direction  vector,  etc.  The  technique  involved 
here  is  a  convergent  one  in  that  after  we  start  up  the  system  by  "assuming" 
a  position  so  we  can  point  the  antennas  and  begin  reading  data,  the  data 
we  obtain  will  yield  position  measurements  that  are  superior  to  our 
initial  assumption.  Once  the  system  is  started  up  it  should  operate 
continuously  even  in  moderately  heavy  seastates.  In  this  way  our 
"assumed"  position  will  be  based  on  past  measurements  and  the  use  of 
the  equations  of  evolution';  and  the  position  should  be  very  accurate. 


(123) 


(124) 
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There  are  six  sources  of  measurement  error  in  this  formulation 
whose  effects  will  have  to  be  considered  in  more  detail  in  a  future 
study  in  order  to  make  a  final  evaluation  of  the  realizable  accuracy 

u 

of  the  interferometer-navigational  system.  These  sources  of  error, 
t  not  necessarily  in  order  of  dominance,  are 

a)  The  effect  of  the  noise  term  nQ(t)  on  the  phase 
measurement  of  the  cosine  term. 

b)  The  accuracy  of  the  ship  motion  measurements  using  the 
accelerometer  technique  described  in  section  2.1. 

c)  The  accuracy  of  the  bearing  angle  measurement,  yQ(t). 

(In  this  worker's  opinion  this  measurement  may  be  the 
dominant  source  of  error  in  the  entire  technique) . 

d)  The  error  induced  by  assuming  a  spherically  symmetric 
stratified  index  of  refraction  in  developing  the 
refraction  matrix  approximation  in  section  2,3  and  in 
developing  the  time-varying  phase  path  difference 

j  formulation  in  section  3.1. 

e)  The  error  in  measuring  the  surface  refractivity  N^Ct,^) 

and  index  of  refraction  n  ( t  -  to  )  used  in  the  refraction 

s  o 

matrix  approximation  and  time-varying  phase  path  difference 
formulation. 

f)  The  accuracy  of  the  ship  velocity  measurement  used  in 

f  ,  f  ,  g  and  g  given  in  equations  (29)  to  (32) . 
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